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ABSTRACT 

We present a study of the vertical magnetic field of the Milky Way towards the Galactic poles, determined 
from observations of Faraday rotation toward more than 1000 polarized extragalactic radio sources at Galactic 
latitudes \b\ > 77°, using the Westerbork Radio Synthesis Telescope and the Australia Telescope Compact 
Array. We find median rotation measures (RMs) of 0.0 ± 0.5 rad m"-^ and +6.3 ± 0.7 rad m~^ toward the north 
and south Galactic poles, respectively, demonstrating that there is no coherent vertical magnetic field in the 
Milky Way at the Sun's position. If this is a global property of the Milky Way's magnetism, then the lack of 
symmetry across the disk rules out pure dipole or quadrupole geometries for the Galactic magnetic field. The 
angular fluctuations in RM seen in our data show no preferred scale within the range w 0?1 to « 25°. The 
observed standard deviation in RM of ~ 9 rad m~^ then implies an upper limit of ~ 1 /xG on the strength of the 
random magnetic field in the warm ionized medium at high Galactic latitudes. 
Subject headings: magnetic fields — ^Faraday rotation — ^polarization — Galaxy: halo 



1. INTRODUCTION 

Large scale coherent magnetic fields are observed in our 
Milky Way and in external galaxies (Beck 2008, 2009); these 
fields play crucial roles in many astrophysical processes in 
the interstellar medium (ISM) — they help to exert pressure 
to balance ordinary matter against gravity and trigger star for- 
mation, they are also responsible for the confinement of cos- 
mic rays, and they can regulate and trace a large scale galac- 
tic wind. Therefore, to better understand galaxy evolution, 
it is necessary to investigate the structure, origin and evolu- 
tion of galactic magnetic fields. The primordial theory and 
the dynamo theory are the two possible explanations for the 
existence of a galactic scale magnetic field. These two the- 
ories make certain predictions on the symmetry of the large 
scale magnetic field with respect to the rotation axis and the 
mid-plane of the galaxy (e.g.. Beck et al. 1996). The charac- 
terization of the overall magnetic field geometry of a galaxy 
allows one to distinguish between various primordial and dy- 
namo models. 

Studies of the large scale magnetic field geometry in the 
Milky Way have been mainly focused on the disk field sym- 
metry with respect to the rotation axis (e.g., Han et al. 2006; 
Brown et al. 2007) even though sight lines towards high 
Galactic latitude may be less turbulent and have less tangled 
fields than are seen in the Galactic plane and hence global 
field patterns can be identified more easily. The strength 
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and the synnmetry of the vertical, azimuthal and radial com- 
ponents of the Galactic magnetic field across the plane can 
provide us with unique insights on the mechanisms which 
maintain the Galactic magnetic field. For example, a rever- 
sal across the Galactic plane, in the azimuthal component of 
the large-scale magnetic field but not in its vertical component 
would indicate a field of dipolar structure, which could re- 
sult from weak differential rotation (near solid body rotation) 
(Ferriere 2005); or a substantial primordial field (Zweibel & 
Heiles 1997). On the other hand, if there is a reversal in 
the vertical but not the azimuthal component across the mid- 
plane, it would indicate a field of quadrupolar structure re- 
sulting from dynamo action due to differential rotation in the 
Galactic disk (Zweibel & Heiles 1997). The dynamo theory 
predicts a weak vertical field compared to the horizontal field 
on galactic scales. If the measured large-scale vertical field 
strength is substantially larger, it could imply a primordial 
component to the Galactic magnetic field (Ruzmaikin et al. 
1988). The ratio of the vertical to horizontal field strength 
also regulates the confinement of cosmic rays to the Galactic 
disk, which can help us better understand the disk-halo inter- 
action in the Milky Way. 

Some information on the large scale structure of the Galac- 
tic magnetic field can come from optical starlight polarization 
and radio synchrotron polarization (see Beck et al. 1996; Beck 
2008, 2009, for a summary). A series of papers: Berdyugin 
et al. (2000); Berdyugin & Teerikorpi (2001); Berdyugin et al. 
(2001, 2004) found that towards flie soufli Galactic pole, the 
polarization of stars traces a field orientation along £ of 80° 
which is parallel to the local spiral arm field, while the same 
field direction could not be traced towards the North Galactic 
pole. Most of our knowledge of the geometry of the large- 
scale Galactic magnetic field to date comes from the Faraday 
rotation measure of distant extragalactic radio sources (EGSs) 
and pulsars. Faraday rotation is a birefringence effect when 
hnearly polarized Ught travels through a magnetized media. 
The plane of the polarization rotates through an angle (in 
radians) given by 

AV; = RMA^ (1) 
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where A is the wavelength of the radiation measured in meters 
and RM is the rotation measure, defined as the integral of the 
Une of sight magnetic field By (in /xG) weighted by the ther- 
mal electron density tieil) (in cm~^) over a line-of-sight line 
element dl (in pc). 

pobserver 

RM = 0.812 / ne{l)Bi\{l)dl rad m'^ (2) 

J source 

The sign of the RM gives the direction of the line of sight 
component of the average field: a negative RM represents 
a magnetic field whose line of sight component is directed 
away from us. With independent knowledge of the thermal 
electron density, one can determine the average Une-of-sight 
magnetic field strength. Faraday rotation is complementary to 
other measurement techniques since RMs provide the direc- 
tion of the magnetic field, while most other techniques (apart 
from Zeeman spUtting) provide the field orientation, but not 
its direction. 

High latitude magnetic fields in the Milky Way have been 
investigated using RMs of EGSs by Morris & Berge (1964), 
Andreasyan (1980), Andreassian & Makarov (1988) and Han 
& Qiao (1994). They all found a strong antisyrmnetric RM 
pattem between the northern and southern Galactic hemi- 
spheres. These authors attribute the pattern to a horizontal 
magnetic field that reverses direction above and below the 
Galactic plane. Han & Qiao (1994) and Han et al. (1999) 
estimated the local vertical magnetic field strength to be ~ 
0.3 /uG, pointing from the south to the north Galactic pole as- 
suming a priori a dipolar field geometry. A wavelet analysis 
on an all sky RM catalogue presented by Frick et al. (2001) 
found that the RM distribution at the largest scales is shifted 
to negative latitudes, indicative of a stronger magnetic field 
in the southern hemisphere and hence the possible existence 
of an antisymmetric^ halo magnetic field^. Unfortunately, lo- 
cal distortions such as that of a Parker instability loop cannot 
be ruled out. A dipolar halo field is also suggested by Sun 
et al. (2008), who demonstrated that an oppositely directed 
torus-Uke halo field above and below the Galactic plane and 
a symmetric disk field provide a reasonable fit to the latitude 
extension of the Canadian Galactic Plane Survey (CGPS) RM 
measurements, but it is unclear if this model can predict RMs 
that are consistent with observations for all ranges of i and b. 

Jansson et al. (2009) found that the halo field might not be 
anti-synmietric globally since a halo field that reverses across 
the plane only towards the inner Galaxy provides a better fit 
to the WMAP 22GHz all sky polarization map and an all- 
sky catalogue of ^ 1,400 RMs of EGSs. Recently, Taylor 
et al. (2009) used EGS RMs calculated from the NRAO VLA 
Sky Survey (NVSS) catalogue, a 1.4GHz continuum survey 
of the entire sky north of DEC=-40°, to derive a vertical field 
strength of 0.30 ± 0.03 /iG for b<0° pointing from the south 
pole towards the north pole and a field strength of 0.14 ± 0.02 
/zG for b>0° pointing from the north towards the south pole. 
This new result is incompatible with the simple dipolar field 
model concluded in earlier studies. 

The current knowledge of the parity of the Galactic mag- 
netic field is based on the sparse all-sky RM measurements 

' In this paper, the terms symmetric and antisymmetric are used to describe 
the vertical and horizontal components of the Galactic magnetic field instead 
of the total magnetic field vector itself. 

' Throughout the paper, we refer to halo fields as non-disk fields that are 
in regions with sufficient diffuse interstellar electrons to produce Faraday ro- 
tation 



and their potential unreliable values because they were de- 
rived using polarization angles at only a few and sometimes 
widely separated wavelengths. Moreover, most conclusions 
were drawn using low-latitude RMs which are likely to be 
contaminated by turbulence and tangled fields in the Galac- 
tic disk. In this paper, we present an accurate rotation mea- 
sure survey of more than 1,000 polarized extragalactic sources 
previously identified in the NVSS catalogue towards both the 
north and the south Galactic poles at a sampling density of 
approximately 1 source/deg'^. Our goal is to measure the 
magnetic field structure at high latitude in great details and 
to study the symmetry of the vertical Galactic magnetic field 
across the disk plane. In § 2, we describe data acquisition 
and reduction procedures. The RM extraction procedure is 
outUned in § 3. The results are presented in § 4. We dis- 
cuss different possible origins of the measured RM pattem in 
§ 5 and derive the corresponding magnetic field properties, 
including both the coherent (§ 5.2) and the random magnetic 
field strength (§ 6). 

2. OBSERVATIONS AND DATA REDUCTION 

To measure the vertical magnetic field at the location of the 
Sun, RMs along many sight lines towards the Galactic poles 
are desired. Since there is only one pulsar with known RM to- 
wards each of the Galactic poles (see § 5.2.1), one has to rely 
on RMs of EGSs to measure the high latitude Galactic mag- 
netic fields. To avoid the inefficiency of conducting a blind 
survey, we applied appropriate filters to the NVSS catalogue 
(Condon et al. 1998) to yield a target list consists of highly po- 
larized compact EGSs near the Galactic poles, towards which 
one can measure RMs. In order to determine RMs as accu- 
rately and precisely as possible, we chose to observe point- 
like extragalactic sources with a signal to noise ratio in linear 
polarization of at least 8. To allow for a reasonable observ- 
ing time per source, we set our polarized flux threshold at 4 
mJy. A rough estimation shows that a vertical magnetic field 
of strength 0.2 /iG and a total column density of thermal elec- 
trons towards the Galactic poles of ^ 25 pc cm~^ (Gaensler 
et al. 2008) would produce a RM of magnitude 4 rad m"^. 
Since the intrinsic RM dispersion of extragalactic sources is 
roughly 15 rad m"^ (Gaensler et al. 2005), one needs approx- 
imately 500 sources to detect the predicted RM signal at a 
decent signal-to-noise level. The polarized flux cutoff at 4 
mJy and the need of ^ 500 EGSs requires us to select sources 
down to a Galactic latitude \b\of 77° in the NVSS catalogue. 

2.1. North Galactic Cap 

RM of EGSs towards the North Galactic pole were acquired 
at the Westerbork Synthesis Radio Telescope (WSRT) over 
the period 2006 July 23rd to August 20th using the 9A=96m 
array configuration spanning baselines from 96 m to 2760 
m; with 8 20-MHz frequency bands (each band is composed 
of 64 channels of width 0.3125 MHz) centered on each of 
1381MHz and 1713MHz. Additional observations were made 
on 2009 April 7th to 8th with the mini-short configuration, 
spanning baselines from 96 m to 2723 m centered on only 
1381MHz. 

We selected 500 EGSs with polarized intensity greater than 
4 mJy at b > +77° from the NVSS catalogue. For each 
source, we obtain 30 cuts of 10 seconds at both 1381MHz 
and 1713MHz (except for sources observed in the additional 
session in 2009), resulting in an integration time of 5 minutes 
per source for each band and hence a total observing time of 
rougly 90 hours. Upon inspection of the data, we found that 
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the IF centered around 1713MHz was heavily contaminated 
by radio frequency interference (RFI) and hence we decided 
to only use the IF centered around 138 IMHz. The overall sen- 
sitivity of each target is ^ 0. 1 mJy-beam"^ . At least one of the 
standard absolute flux calibrators 3C286, 3C147 or CTD93 
was observed before and after each 12 hour observing run 
to allow for amplitude and phase calibration. Bandpass and 
leakage correction was determined by observing an unpolar- 
ized source, either 3C147 or CTD93. A source with known 
polarization, one or both of 3C286 and 3C138 was observed 
for absolute polarization angle calibration. 

The NRAO's Astronomical Image Processing System 
(AIPS) was used for data edition and calibration. AIPS is 
designed to process interferometric data for telescopes with 
circular feeds while the WSRT consists of orthogonal Unear 
feeds (Weiler 1973). Therefore, to properly calibrate and re- 
duce WSRT data in AIPS requires various non-standard pro- 
cedures, such as relabeling the Stokes parameters. These 
caUbration steps are described in details in the AIPS reduc- 
tion cookbook for WSRT data^. After calibration, flagging 
and rebinning the 512 0. 3125MHz- wide channels results in 
16 10-MHz wide channels centered around 20 cm. The u- 
V data were then exported to FITS format so that imaging 
could be carried out in the MIRIAD package (Sault & Killeen 
2003). For each pointing and for each of the 16 10-MHz 
wide frequency channel, maps of Stokes parameters Q and 
U were made using the INVERT task in MIRIAD with natu- 
ral weighting to maximize sensitivity. These maps were then 
deconvolved using the CLEAN algorithm with a threshold of 
5 fT. A restored map for each pointing, for each Stokes Q and 
U at each of the 16 channels around 20 cm were made. This 
results in a total of ^ 16,000 channel maps, each with a sen- 
sitivity of roughly 0.4 mJy-beam"' and resolution of ^ 0.5'. 
A linearly polarized intensity (PI) map, corrected for positive 
bias was made for each source. The brightest polarized pixel 
of each EGS is identified and its Stokes Q and U values across 
the frequency band were extracted for RM determination. 

2.2. South Galactic Cap 

RMs of EGSs towards the south Galactic pole were ac- 
quired at the Australia Telescope Compact Array (ATCA) 
over the period 2006 June 23rd to July 4th, using the 1 .5D 
array configuration spanning baselines from 107.1 m to 1439 
m (without antenna 6), and 2008 June 22nd and June 25th, 
using the EW352 array configuration spanning baselines from 
30.6 m to 352 m (without antenna 6). All ATCA data consist 
of 32 adjacent frequency channels each of bandwidth 4 MHz 
centered on each of 1384 MHz and 2368 MHz. 

The standard primary flux and bandpass caUbrator PKS 
B 1934-638, whose flux at 1384 (2368) MHz was assumed 
to be 14.94 (11.14) Jy, was observed at the beginning and 
the end of each observation. A secondary calibrator (PKS 
2339-353, PKS 0023-26, PKS 0010-401 or PKS 0153-410) 
was observed every hour and was used to correct for polar- 
ization leakages and to calibrate the time-dependent antenna 
gains (Sault et al. 1991). We have applied the same EGS se- 
lection criteria as that towards the north pole and selected 533 
sources towards the south Galactic pole. Allowing ~ 15 min- 
utes per pointing results in a total observing time of approxi- 
mately 150 hours. The overall sensitivity of each field is ^ 0.3 
mJy-beam~' . Using the same pipeline that processes the north 

' http://www.astron.nl/radio-observatory/astronomers/analysis-wsrt- 
data/analysis-wsrt-dzb-data-classic-aips/analysis-wsrt-d 



Galactic pole data, we have imaged EGSs towards the south 
Galactic pole, resulting in images of resolution ^ 0.5', simi- 
lar to that of the north cap data. A total of ~ 26,500 Stokes Q 
and U channel maps (13 channels around 13 cm and 12 chan- 
nels around 20 cm) were produced, each with a sensitivity of 
roughly 1 mJy-beam"'. 

2.3. Ha Observations 

The warm ionized medium (WIM), consisting of diffuse 
ionized gas, is a major constituent of the interstellar medium 
in the Milky Way. The WIM can be probed using the Ha 
line and other optical recombination lines. The Wisconsin 
Hq Mapper (WHAM) is a 1° -angular resolution, velocity re- 
solved survey of the Ha line in the northern (completed) and 
the southern (in progress) sky (Haffner et al. 1999). Ha mea- 
surements together with dispersion measure (DM) 

DM = / Hgdl = EgL pc cm"^ (3) 
Jo 

of pulsars along a line of sight can provide information on 
the dumpiness of the WIM (see for example, Reynolds 1991; 
Gaensler et al. 2008). The observed Ha intensity in Rayleighs 
(1 R = 10'' photons per 4tt steradian = 2.42 x 10"^ ergs cm"^ 
s~' sr~') towards the Galactic caps can be converted into emis- 
sion measure (EM) assuming a constant electron temperature 
around 10,000K and negUgible extinction at high Galactic lat- 
itudes (< 0.02 mag Schlegel et al. 1998) using the following 
equation (Reynolds 1991) 

EM = 2.75(7/ 10Y%a. (4) 
The EM along a sight Une is defined as 

EM= / n^Jl = r^^L pccm-^, (5) 
Jo 

where is the average of the electron density squared along 
the total path length L. 

The WHAM survey has fuU coverage of the north Galactic 
cap (b>+n°) and currently a ~ 50% coverage of the south 
Galactic cap (b< -77°). The average surface brightness is 
about 0.5 R(EM ~ 1 pc cm"^) towards the caps. We have 
verified that there are no discrete HII regions at |fo|>77°. The 
WHAM data towards the cap regions are used to search for 
correlations between EM and RM in § 4.1. 

3. RM COMPUTATION 

Traditionally, the RM of an EGS is calculated by a linear 
least square fit of the unwrapped polarization position angle 

ih= - tan"' — (6) 
2 2 

as a function of the square of the wavelength. However, since 
the polarization angle can wrap around nn, there exists mul- 
tiple RMs that provide equally good least square fits to the 
observed angle versus relation. We have used the Faraday 
rotation measure synthesis, first presented by Bum (1966) and 
more recently illustrated in Brentjens & de Bruyn (2005) and 
Heald et al. (2009), to obtain reliable RMs of the observed 
EGSs. This method is equivalent to taking the Fourier trans- 
form of the complex polarization vector measured at each fre- 
quency channel. The RM Synthesis method does not operate 
with polarization position angles and hence does not lead to 
the rnr ambiguity problem. Also, non-Gaussian statistics of 
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polarization angles can be avoided. Moreover, this method 
can resolve multiple RM components along the line of sight 
or within the telescope beam, for which case the angle ver- 
sus fit would not give a sensible result. Backends of the 
ATCA and the WSRT make the implementation of RM syn- 
thesis possible as Stokes Q and U information can be obtained 
in numerous neighboring frequency channels. 

We have performed RM Synthesis on the observed EGSs 
following the method described in detail by Brentjens & de 
Bruyn (2005). Similar to an aperture synthesis experiment 
where large gaps in the uv-plane causes high-level sidelobes 
in the image plane, incomplete wavelength coverage of the 
RM experiment leads to sidelobes in the resulting Faraday 
depth spectra. 

This can cause serious problems when attempting to ex- 
tract the correct RM. To enable accurate extraction of both 
the RM and the polarized flux, we have deconvolved the com- 
plex Faraday depth spectrum using the CLEAN algorithm 
as described by Brentjens (2007) (also see the Appendix in 
Heald et al. 2009). We have adopted a gain factor of 0.1 
and we stop CLEANing once the residual peak falls below 
4 times the noise level in the final spectrum. The value of 
RM and the polarized intensity of the EGS is calculated by 
fitting a parabola to the peak(s) of the deconvolved Faraday 
depth power spectrum. The precision of a RM measurement 
is determined by the signal-to-noise of the polarization detec- 
tion and the FWHM of the rotation measure transfer function 
(Equation (61) in Brentjens & de Bruyn (2005)). The mea- 
surement uncertainty of an RM with a debiased peak polar- 
ized intensity'' in the Faraday depth spectrum of PIpeak. a rms 
error in the CLEANed Faraday depth spectrum a and a dirty 
beam FWHM of 5(j) rad m~^ is 



ARM: 



2PIpeak/ 0" 



(7) 



The FWHM of the dirty beam is 343 rad m'^ for the WSRT 
observations while the value for the ATCA observations is 99 
rad m"^. The ATCA observation has relatively low sensitivity 
but small 64), while the WSRT observation has better sensi- 
tivity but larger 54). This results in comparable RM errors for 
both the northern and southern data sets. We have verified, 
using several sources with steep/ shallow spectra in our sam- 
ple, that the spectral index has little effect on the location of 
the Faraday dispersion function peak (Brentjens & de Bruyn 
2005). We therefore choose not to include spectral index ef- 
fect in our RM computation/ deconvolution procedure. 

We have followed the above steps to derive RMs of EGSs 
towards the north Galactic cap using the 20 cm WSRT data. 
For sight lines towards the south Galactic cap, RMs were first 
computed using only the 20 cm ATCA data. If the source 
appears to have a well behaved Faraday depth spectrum with 
a single RM component, we then add in the 13 cm ATCA 
data to achieve better RM precision. This step requires the 
assumption that the same RM component dominates at both 
13 and 20 cm. To minimize the potential systematics (such 
as caUbration errors) introduced when we calculate the RMs 
using both bands, we have rejected sources that have RM val- 
ues derived from the 20 cm only data that are at 2a level dif- 
ferent from those obtained using both bands. We have also 
discarded sources with polarization detection signal-to-noise 

^ The polarized intensity is debiased to the first order by subtracting the 
noise in the Faraday depth spectrum from the measured peak polarized inten- 
sity in quadrature (Simmons & Stewart 1985). 



ratio below 8 and with polarized fraction greater than 30 
The least square fit and Faraday depth spectra plots of two ex- 
ample EGSs are shown in Fig 1. Since the frequency setup 
for the ATCA observations leaves a larger gap in space 
than the WSRT observations, it leads to higher side lobe lev- 
els seen in the ATCA Faraday depth spectrum. For bright po- 
larized sources (as the two example EGSs), the least square fit 
and the RM synthesis method should yield converging results. 

As the Faraday rotation phenomenon is a line-of-sight ef- 
fect, the observed RM is the sum of aU the RM contributions 
along a sight line. This includes the intrinsic RM which orig- 
inates in the magneto-ionic environment at the radio source; 
the RM produced in the intergalactic medium (IGM) or any 
intervener(s); that produced by the Milky Way's large scale 
and small scale field; as well as that produced by the iono- 
sphere. The RM contribution from the IGM is negligible - 
cosmological Faraday rotation measure has an estimated up- 
per limit of ^ 2 rad m"^ assuming that the cosmological mag- 
netic field and electron density is homogenous (Vallee 1990). 
On the other hand, if these two quantities vary on scales much 
smaller than the Hubble scale, then no cosmological imprint 
to the all sky RM is expected as it would average out to be 
zero. 

As Global Positioning System (GPS) data were not avail- 
able to monitor the real time ionospheric activity during the 
observations at the ATCA and the WSRT, we can only rely on 
theoretical models to estimate the contribution of ionospheric 
RM to the observed RMs. We note that theoretical predictions 
do not take into account of sudden and temporal ionospheric 
changes. We have used the Intemational Geomagnetic Refer- 
ence Field (IGRF) Model^ (MacMillan & Maus 2005), along 
with the International Reference Ionosphere (IRl) 2007^ (Bil- 
itza & Reinisch 2008) to estimate the ionospheric RM during 
the observations. It has been found that the ionospheric RM 
above the ATCA is negative (geomagnetic field points away 
from the center of the earth) and that above the WSRT is posi- 
tive (geomagnetic field points towards the center of the earth) 
with a magnitude < 1 rad m~^, which is consistent with the 
estimation reported by Tinbergen (1996). The variation of 
ionospheric RM along a certain RA and DEC over the course 
of the observation is also < 1 rad m~^. This can be verified 
by computing the RM of a bright polarized EGS as a function 
of time over the observing run. We find that the RM stays 
constant as a function of time within its error To account for 
the RM uncertainties introduced by Faraday rotation from the 
ionosphere, we have added a systematic error of 1 rad m~^ in 
quadrature to the RM measurement error for all sources. 

4. RESULTS 

The RM derivation procedures described in the previous 
section produced 472 reliable RMs towards the north Galactic 
pole (94 % of the total number of observed sources) and 341 
(72 % of the total number of observed sources) towards the 
south. The source coordinates, RMs, RM uncertainties, and 
their flux information are listed in Table 1 and 2 for sources 
towards the north and the south cap respectively. To help vi- 
sualizing the distribution of RMs towards the Galactic caps, 
we present the histograms of the unfiltered RMs in Fig 2. 
They are shown as solid outlined histograms which appear 



Sources with unusually high polarized fraction might be of Galactic ori- 



gin 



* Available at http://omniweb.gsfc.nasa.gov/vitmo/igrf_vitmo.html 
^ Available at http://omiiiweb.gsfc.nasa.gov/vitmo/iri_vitmo.html 
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to be Gaussian with few outliers. The spatial RM distribu- 
tion towards the north and south Galactic caps are presented 
in the top and bottom panel of Fig 3 respectively. In this sec- 
tion, we will address how outliers and anomalous RM regions 
were identified. We do so to ensure that the RMs towards the 
Galactic caps are representative of the large scale Milky Way 
foreground with minimal contributions from elsewhere along 
the sight line. 

4.1. Features towards the Galactic poles at other 
wavelengths 

We have visually inspected multi-wavelength images to- 
wards the Galactic caps to look for possible features seen at 
other wavelengths that might be associated with RM struc- 
tures. In particular, we have examined the following maps: 
the HI neutral hydrogen map (Kalberla et al. 2005); the dust 
map (Schlegel et al. 1998); the ROSAT all sky survey X- 
ray map (Snowden et al. 1997); the 408 MHz total syn- 
chrotron intensity map (Haslam et al. 1982); the 1.4 GHz 
polarized synchrotron emission map (north: WoUeben et al. 
2006, south: Testori et al. 2008), and the Ha all sky compos- 
ite map (Finkbeiner 2003). 

We note that the Coma cluster (the brightest X-ray feature 
in the ROSAT all sky X-ray image towards the north Galactic 
cap) is located only 2° away from the Galactic pole. Accord- 
ing to Kim et al. (1990), the Coma cluster has a magnetic field 
strength of 2 /iG and it can produce an excess RM of sev- 
eral 10s of rad m"^. Hence, we have excluded 10 EGSs within 
2° of the Coma cluster (RA=12:59:48.7, DEC=H-27:58:50.0) 
to ensure that our RM sample is free of contributions from 
the intragalactic medium of Coma. These rejected EGSs are 
denoted by superscript a in column 8 of Table 1 . 

Towards the north Galactic pole, part of the North Polar 
Spur (MPS), a bright polarized feature at radio wavelengths 
and a prominent X-ray feature at high positive Galactic lati- 
tude, is in our field of view. As pointed out by Frick et al. 
(2001), the NFS can produce a maximum foreground RM of 
^ -28 rad m"^ towards £=289°, b=+2'i° , which corresponds 
to an internal magnetic field strength of ~ 0.9 /iG (Frick et al. 
2001; Heiles et al. 1980). This internal magnetic field strength 
of the NFS coupled by thermal electrons can produce sig- 
nificant RMs in our field of view . The polarized emission 
of the NFS and its counterpart below the Galactic plane in 
the DRAO 1 .4 GHz polarization survey has been modeled by 
WoUeben (2007) with two overlapping spherical synchrotron 
emitting magnetized shells where the Sun resides in one of 
them. The model correctly predicts the morphology of the 
NFS polarized emission with the projection of the overlapped 
region onto the sky roughly coincides with the NFS X-ray 
emission. Although the two-shell model does not include pre- 
dictions of RMs through the shells, one expect complicated 
RM structure for sight lines that pass through both shells. 
Therefore, we choose to discard all EGSs whose sight line 
penetrate the polarized radio feature and the prominent X-ray 
feature of the NFS (^<h-80° and 50° <£<3 10°). A total of 56 
EGSs are rejected and they are indicated by superscript b in 
column 8 of Table 1 . At high negative Galactic latitude, there 
are no features seen at other wavelengths that might contami- 
nate RMs. 

One might expect some correlation between RM and Ha as 
they both depend on the electron density content along a sight 
line. Fig 3 shows the distribution of RMs towards the Galactic 
cap regions overlaid on the Finkbeiner (2003) Ha all sky com- 
posite map. No evidence of correlations between RM and EM 



could be seen. We have further investigated this issue by mak- 
ing a scatter plot of EM versus RM towards the north and the 
south Galactic cap (Fig 4). One can clearly see that points in 
these plots appear to be randomly located without any trend. 
The linear correlation coefficient between RM and VEM is 
~ 0.04 towards the north Galactic cap and ^ 0.03 towards 
the south Galactic cap hence there is little support for corre- 
lations between EM and RM towards the Galactic poles. It is 
not surprising, however, to find RM and EM uncorrelated. As 
EM is proportional to the electron density squared, it is more 
sensitive to high density regions; while RM is proportional to 
just the electron density, thus it has more contributions from 
diffuse low density regions that occupy greater Une-of-sight 
path lengths. Magnetic field reversals along the line of sight 
can also remove any possible correlations between RM and 
EM. All these effects work towards decoupling the measured 
RM and EM along a particular sight line. 

4.2. Outliers 

Besides excluding features at other wavelengths that might 
contribute to the observed RMs, one also has to remove 
sources in the sample with high intrinsic RMs. Typical radio 
galaxies at a few GHz are generally dominated by their bright 
radio lobes which are thought to have small intrinsic Fara- 
day rotation measures (see for example, Willis & Strom 1978; 
Feain et al. 2009). Since we have selected EGSs based only 
on their polarized flux and \b\, our sample inevitably contains 
core-dominated EGSs which have high internal RMs (possi- 
bly generated from a cocoon filled with magnetic fields and 
thermal electrons) and EGSs which lie along line of sights to 
multiple intervening galaxies or clusters with high RMs by 
chance. We have employed the Chauvenet's criterion^. For a 
sample size of ^ 400, one expects less than 0.5 source with 
an RM that deviates more than 3.2 a from the median. EGSs 
with RM values 3.2 ct away from the median of the distribu- 
tion are discarded^. This is based on the assumption that the 
underlying RM distribution is Gaussian. The thick vertical 
lines in Fig 2 indicate boundaries beyond which RMs have 
been rejected. These EGSs are indicated in column 8 of Ta- 
ble 1 and 2. 

4.3. Smooth RM maps and anomalous RM regions 

After discarding EGSs in the cap regions which have large 
intrinsic RMs, subregions towards the Galactic caps where 
RMs differ systematically from the overall smooth behav- 
ior expected from a large scale magnetic field should also 
be identified. To locate these subregions, we have first con- 
structed smooth RM maps towards the Galactic caps using the 
region-filtered (without the NFS and Coma Cluster regions) 
and outlier-free RM sample. We have divided the cap regions 
into 2° X 2° cells and computed the median RM of the EGSs 
within the individual cells. To ensure the statistics are reliable, 
we require at least 3 or more EGSs in a cell to derive a valid 
median RM. The resulting smoothed maps are shown in Fig 5. 
One can see clearly from the smoothed maps that the RM to- 
wards the north Galactic pole is distributed randomly around 
zero whereas the overall RM towards the south Galactic pole 

^ The Chauvenet's criterion states that a data point can be discarded with 
reasonable confidence if less than half an event is expected to be farther from 
the median than the suspect point (Bevington & Robinson 2003). 

' The Chauvenet's criterion is applied to the RM data set towards the north 
cap after regions coinciding with the Coma cluster and the NFS have been 
removed. 
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is positive. Next, we have calculated how many standard de- 
viations the median of RMs in a cell deviates from the me- 
dian of the entire RM sample. We discard all sources within 
a cell if the cell has a minimum of 3 EGSs and has a median 
RM that deviates more than 1.65 a (> 90%) away from the 
overall median of the RM distribution. We have identify 1 
such cell towards the north Galactic cap and it is indicated by 
green rectangle in the top panel of Fig 5. No anomalous RM 
cell was identified towards the south Galactic cap as shown 
in the bottom panel of Fig 5. We are aware that this rejec- 
tion criterion is only sensitive to anomalous RM regions with 
scales comparable to the 2° x 2° cells and that the boundary 
of the cells is chosen arbitrarily. We have verified that sources 
being rejected by this algorithm are very similar regardless of 
the choice of the cell size. This leads us to beUeve that the 
method outlined above is a robust way of rejecting anoma- 
lous RM regions. The 3 sources towards the north Galactic 
cap discarded by applying the above criterion are indicated in 
colunm 8 of Table 1 . After the rejection procedures described 
above, there remains 400 RMs towards the north Galactic cap 
and 329 RMs towards the south Galactic cap. These RMs 
should well represent the large scale magnetic field towards 
the Galactic poles. 

4.4. Rotation Measure structure function 

To study the RM structures towards the Galactic caps and 
to understand the properties of turbulence at high |fc|, we have 
computed the second-order ID RM structure function (SF) 
defined as 

5FRM,obs('-) = {{RM{x)-RM{x+ r)f) , (8) 

where r is the projected angular separation between a pair of 
extragalactic sources and angular brackets denote the expec- 
tation value which involves taking the ensemble average of 
independent measurements with the same range in angular 
separation r. The uncertainties of RM measurements (ctrm) 
contribute to the observed RM structure function in the form 
of a DC offset (see Appendix A of Haverkorn et al. 2004). 
The offset-corrected structure function is given by 

SF^ir) = 5FRM,obs('-) -SF^^ir). (9) 

Since RMs of different components of the same EGS will 
contribute many square differences but they do not provide 
independent information (Minter & Spangler 1996), we have 
identified EGSs in the filtered RM data sets that appear to 
be associated with each other using the NASA/IPAC Extra- 
galactic Database (NED) and excluded them from the struc- 
ture function calculations. We defer the study of RM fluctu- 
ation of these closely spaced EGSs to a forthcoming paper. 
The RM structure functions were constructed using 354 RMs 
towards the north and 319 RMs towards the south Galactic 
cap. This is the largest and most reliable RM data set avail- 
able to date to study turbulence at high Galactic latitudes - the 
high source density allows us to probe a large range of scales 
from ~ 0. 1 ° up to 26 ° , while the large number of RMs allows 
one to determine the SF accurately based on 62,481 correla- 
tions towards the north and 50,721 towards the south Galactic 
poles. 

The resulting RM SF towards the north and the south Galac- 
tic caps are plotted in the top and bottom panel of Fig 6 respec- 
tively. A power-law fit is performed on the observed RM SFs 
- the slope of the SF towards the north Galactic cap is +0.08 
± 0.01 while that towards the south is +0.03 ± 0.01. Since 
the structure function measures the fluctuation of RM along 



different sight lines, it is not sensitive to constant RM contri- 
butions in the cap regions, such as that from a constant vertical 
magnetic field. However, the RM SF is sensitive to any large 
scale magnetic field or electron density gradients in the field 
of view. For example, a large scale horizontal magnetic field 
towards the Galactic poles would produce more structures on 
the largest scales and would lead to a rising structure function. 
The nearly flat structure functions towards both caps confirm 
our visual inspection of Fig 5 that RMs are uncorrelated and 
that there is no preferred scale of RM structure towards the 
Galactic caps. The flat RM SFs towards the Galactic poles 
are consistent with the non-detection of horizontal fields to- 
wards the poles in § 5.2.2. 

4.5. The RM distribution towards the Galactic poles 

The distribution of the final RM data sets is shown in the 
shaded histograms in Fig 2. The median RM towards the 
north Galactic pole is 0.0 ± 0.5 rad m~^, which indicates that 
the median RM towards the north Galactic cap is consistent 
with zero. The 3 a upper limit on the |RM| towards the north 
Galactic pole is 1.4 rad m~^. On the other hand, the median 
of RMs towards the south Galactic cap is +6.3 ± 0.7 rad m"^. 
The RM distribution towards the south Galactic cap is incon- 
sistent with zero at greater than 9 a level. 

4.6. Comparison with Taylor et al. (2009) 

As mentioned in § 1, Taylor et al. (2009) have computed 
RMs from the NVSS catalogue. In Figure 7, we have plotted 
the RMs derived by Taylor et al. (2009) against the RMs that 
we have derived for the same sources. The thick solid line of 
slope 1 in the figure indicates where sources should lie if the 
NVSS RMs and our RMs are equal. Approximately 57% of 
the RMs from the two samples towards the north cap agree 
with each other within their measurement errors, whereas this 
percentage is 53% towards the south cap'". We do not find 
the ± 650 rad m"^ RM ambiguity problem as discussed by 
Taylor et al. (2009) in our high latitude data set because the 
measured RMs have small magnitudes. The linear correlation 
coefficients between the two RM data sets are: 0.39 towards 
the north and 0.36 towards the south, indicating that there ex- 
ists substantial difference between the two RM samples. The 
signs of the Taylor et al. (2009) RMs and our RMs do not al- 
ways agree - this is demonstrated by sources located in the up- 
per left and lower right corners of the scatter plots. The differ- 
ence between Taylor et al. (2009) RMs and our WSRT/ ATCA 
RMs is most likely due to multiple RM components that pro- 
duces non-linearity in polarization angle against relation 
that Taylor et al. (2009) was not able to identify, and poten- 
tially ionospheric RM which depends on the observing time, 
site and also epoch in the solar cycle. While RMs derived in 
Taylor et al. (2009) can be used collectively to describe the 
large scale Galactic magnetic field by averaging over large ar- 
eas, extra care must be taken if one plans to use the individual 
RM values in their catalogue since these RMs can potentially 
be inaccurate. Our RM catalogue is more suitable if one plans 
to investigate the RM properties of individual sources. 

5. THE ORIGIN OF THE OBSERVED HIGH GALACTIC LATITUDE 
RM STRUCTURE 

In this section, we explore possible origins of the observed 
RM pattern towards the Galactic caps: an overall RM consis- 

One would expect approximately 68% of the RMs from the two studies 
to be consistent within I it if there were no systematic errors. 
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tent with zero towards the north and a positive median RM to- 
wards the south. We compare predictions from different mod- 
els with the observed RM pattern. We first investigate local 
sources/ events that might give rise to the observed RMs: the 
local interstellar medium, the local bubble and a Parker insta- 
bility loop". We then estimate the magnetic field strength to- 
wards the Galactic poles from the observed RMs by assuming 
that the observed Faraday rotation occurs in diffuse ionized 
gas in the Galaxy. Finally, we consider the likelihood that the 
observed RM pattern has been generated by global events in 
the Galaxy, such as a large scale Galactic wind, large scale 
Galactic dynamos or a relic field. 

5.1. Local origins 

5.1.1. Local Interstellar medium 

The immediate medium surrounding the sun could poten- 
tially produce the observed RM pattern at high Galactic lat- 
itude. Even though the sun resides in a low density cavity 
called the Local Bubble (LB) (see more in § 5.1.2 ), there are 
warm partially ionized clouds surrounding the Sun (Linsky & 
Redfield 2007). These clouds are magnetized and have free 
electrons, and therefore are capable of rotating the polariza- 
tion plane of incident radiation. Spangler (2009) estimated the 
upper limit of |RM| produced by the LISM to be 0.32-1.1 rad 
m"^ by assuming that typical clouds have an electron density 
of 0.12 cm-3, a volume filling factor of 5.5% - 19% (Red- 
field & Falcon 2008; Redfield & Linsky 2008) and a mag- 
netic field strength of 4 /xG. Different studies have yield dif- 
ferent estimates for the local magnetic field strength. Snow- 
den et al. (1998) suggested that a magnetic field of up to 7/iG 
is required in the LB to counter balance the enormous ther- 
mal pressure {pjk^ 15,000 cm~^ K ) exerted by the enclosing 
hot X-ray gas, but the recent discovery of X-ray emission as- 
sociated with charge exchange between solar wind ions and 
heliospheric plasma has greatly alleviated the need of non- 
thermal pressure support of the LB and has lowered the re- 
quired magnetic field strength to ^ 2.8 /xG (Welsh & Shelton 
2009). Other works (Opher et al. 2007; Wood et al. 2007) 
find a local field strength of ~ 2 /zG. Since various estima- 
tions of the local magnetic field strength are lower than the 
4 ^G adopted by Spangler (2009), the LISM |RM| can be a 
factor of two smaller than the estimates of Spangler (2009). 
We note that the estimated LISM |RM| contribution of 1 
rad m"^ is an upper limit as |RM| would be much smaller if 
the magnetic field reverses within the clouds or from cloud to 
cloud. As the typical RM measurement error in our experi- 
ment is a few rad m"^, contribution of RM from the LISM is 
likely to be negligible. 

5.1.2. Effects of the local bubble wall 

In this section, we consider the possibiUty that the Faraday 
rotation originates in the wall of the Local Bubble. As men- 
tioned in § 5.1.1, the Sun is situated in a low-density cavity 
thought to be created by star formation and subsequent super- 
nova explosions that occurred in the past 25-60 Myrs (Frisch 
2007). Such processes would sweep up materials and mag- 
netic fields in the solar neighborhood into a dense shell and 
might produce measurable RM for sight lines through it. 

To test if the LB wall is responsible for producing the ob- 
served RM towards the Galactic caps, we have adopted the 

' ' We have excluded the NPS two shell model (Wolleben et al. 2006) from 
consideration as we have already discarded RMs of EGSs whose sight lines 
intercept the NPS (see § 4.1) 



LB wall model constructed by Cordes & Lazio (2002, 2003) 
with slight modifications: we have modeled the LB wall as 
a slanted cylinder of constant radius 0.085 kpc that extends 
to 0.2 kpc both above and below the Galactic plane. Unlike 
the original model in Cordes & Lazio (2002, 2003), our cylin- 
drical wall has open ends since the recent expansion of Nal 
measurements made by Lallement et al. (2003) and Welsh 
et al. (2009) find no continuous neutral LB boundary at high 
Galactic latitudes. We have computed the projection of the 
modeled LB boundary towards the Galactic poles and plotted 
it in Fig 3 as a green dotted line. None of the sight lines to- 
wards the north Galactic cap intercepts the local bubble wall, 
while 267 out of 341 sight lines intercept the modeled wall 
towards the south Galactic cap. Towards the south Galactic 
pole, the median RM of sight fines that penetrate the LB wall 
is -1-6.1 ± 0.7 rad m~^ whereas those do not penetrate the LB 
wall have a median of + 7.8 ± 1.8 rad m"^. Since the inferred 
RM through the LB wall is consistent with zero (-2 ± 2 rad 
m"^), we conclude that the local bubble wall is unlikely to be 
the major contributor to the observed RM towards the caps if 
the adopted wall model is realistic. 

The 3(7 upper limit of |RM| ^ 7 rad m"^ through the local 
bubble wall can be used to infer the magnetic field strength 
in the LB wall. Since Bhat et al. (1998) found that the scin- 
tillation measures of 20 nearby pulsars were well modeled 
by a scattering structure with local electron density enhance- 
ment of a factor of 10 which roughly coincides with the neu- 
tral LB wall, we assume that the electron density wall can be 
well traced by the neutral wall. If energetic events have swept 
up magnetic fields and electrons of initial density 0.025 cm"^ 
(Bhat et al. 1998) into a cylindrical shell of radius /? ~ 85 pc 
(Cordes & Lazio 2002, 2003) centered at the Sun, the thick- 
ness of the shell in pc {5R) is related to the electron density 
enhancement (x=10) in the shell and its radius R by conserv- 
ing the mass of electrons before and after the formation of the 
shell hy 5R = j^. The thickness of the cylindrical LB shell 
is estimated to be ^ 4 pc. If the RM produced by the wall 
is 7 rad m"^, it implies a magnetic field strength of ^ 9 /xG 
. This prediction from a simple theoretical model of the LB 
is in rough agreement with the magnetic field strength of 8 
HG estimated by Andersson & Potter (2006) (or earUer Leroy 
1999), who have applied the Chandrasekhar & Fermi (1953) 
method on starlight polarization measurements towards stars 
at distances from 40-200 pc in the direction of £=300°, b=0° . 

5.1.3. Small Scale outflows: a Parker instability loop 

Parker (1979) demonstrated that a system of horizontal 
magnetic field and cosmic rays in a vertical gravitational field 
can be unstable with respect to the bending of magnetic field 
lines in the vertical plane. For example, energetic stellar 
events can produce such waviness in a Galactic disk (Kron- 
berg 1994). It is possible that the entire surface of the Galac- 
tic disk is packed with Parker loops with height ^ 1 kpc and 
width of 0.1-1 kpc (Parker 1992). If these magnetic loops 
thread the warm ionized gas in the Galaxy, it is plausible 
that they contribute towards the observed RM patterns at high 
Galactic latitude. 

Frick et al. (2001) performed a wavelet analysis on the all 
sky RM distribution and found that the RM structure at the 
largest scale has been shifted to a negative Galactic latitude 
of -15°, even after omitting sight lines that intercept Loop I. 
While one can interpret this shift as due to a stronger large- 
scale field in the southern Galactic hemisphere, possibly due 
to a separate halo dynamo with opposite parity to the disk 
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field, it can also be due to a Parker instability loop. The au- 
thors suggest a scenario that a Parker instability loop with the 
Sun located near its top (~ 50 pc above the galactic plane), 
a horizontal extent of ~ 400 pc and a magnetic field strength 
enhancement of 0.5 fj,G can shift the symmetry axis of the 
RM structure to b=—l5°. While the physical parameters of 
the inferred loop is consistent with that proposed by Parker 
(1992), it has difficulties explaining the observed RM pattern 
towards the Galactic caps on its own. Since the sun is located 
near the top of this loop, when looking towards a cone of ra- 
dius 13° around the Galactic poles, the extra |RM| produced 
by the loop is small ^1-2 rad m"^, comparable to individ- 
ual RM measurement errors. Furthermore, the median RM 
towards the south Galactic cap would be zero if the loop is 
sufficiently symmetric as half of the south cap region should 
have the exact opposite RMs to the other half. This is in- 
consistent with what has been observed: the RM towards the 
entire south Galactic cap is positive (Fig 5). If the Sun is not 
located exactly at the top of this loop, then one would expect 
median RMs of different signs towards each cap because the 
magnetic field lines should be continuous. This is again in- 
consistent with the observations because the average RM to- 
wards the north cap is consistent with zero. We therefore con- 
clude that the Parker instability loop proposed by Frick et al. 
(2001) alone cannot reproduce our high latitude RM measure- 
ments. 

We note that the argument above is based on the assumption 
that Parker instability loops exist in the warm ionized gas (typ- 
ical WIM density ~ 0.1 cm"^) in the Galaxy. If these instabil- 
ities are associated with the hot phase of the ISM where the 
typical magnetic field strength is ~ 0.1 //G (Beck et al. 1996) 
and the typical density is very low 10""' - 10"^ cm"-' (Sembach 
et al. 2003), the expected RM from such loops would be much 
smaller than the observed |RM| of a few rad m"-^. In this case, 
such loops could not produce the RMs seen towards the caps. 

5.2. Large scale magnetic field in the halo of the Milky Way 

In this section, we consider the possibility that the observed 
RM pattern originates on the Galactic scale. We first esti- 
mate the implied magnetic field strength and direction from 
the observed RMs assuming that the observed Faraday rota- 
tion occurs in diffuse ionized gas. In reality, it is mostly in the 
WIM that the Faraday rotation takes place because even with 
a higher filling factor, the hot halo electrons (~ 10^ - 10^ K) 
are of very low densities (< 10"^ - 10"^ cm"-'). Therefore, it 
is reasonable to assume that the hot halo electrons have neg- 
ligible contributions to the observed RMs. We then consider 
the possibiUties that the observed magnetic field is the result 
of either a Galactic wind, a large scale Galactic dynamo or a 
primordial field. 

5.2. 1 . Vertical Magnetic Field 

After the removal of extreme RMs and anomalous RM re- 
gions from the data set following the method described in § 4, 
one can infer the properties of the Galactic magnetic field 
from the remaining EGSs RMs assuming that the WIM is 
where most of the Faraday rotation takes place. 

The observed RM is the integral of the projection of the 
Galactic magnetic field along the fine of sight weighted by 
thermal electron density. At the highest latitude, RMs mea- 
sure mostly the vertical component of the Galactic magnetic 
field, as the projection of the horizontal component along the 
line of sight is very small. In a right handed coordinate system 



centered at the location of the Sun where the positive z axis 
points from the south Galactic pole to the north Galactic pole 
and the positive x axis points towards the Galactic center, the 
Galactic magnetic field is of the form 

B = Bh cos{k)x+BH sm{lQ)y+B^z, (10) 

where B^ is the coherent component perpendicular to the 
Galactic disk (defined such that a positive B^ implies a field 
pointing from the south to the north Galactic pole), and Bh is 
the coherent component parallel to the Galactic disk directed 
along Galactic longitude Iq. We consider B^, Iq and Bh to 
be constant within the 13° radius cone that we are probing 
around the Galactic poles. If we use the definition of RM 
(Eq 2) and assume that the volume averaged thermal electron 
density in the Milky Way is an exponential disk of mid- 
plane density nefi and scale height Hq : 

«,(z) = n,,oe"l^l/^», (11) 

then the RM towards Galactic coordinates {t,b) can be ex- 
pressed as 

RM = QMlnefiHQ{aB^-BHCOs{t-to)/X3isv\b\), (12) 

where a =+l for b<-ll° and a= -1 for b>+ll° . 

If we take the average RM, (RM), along many sight lines 
with full coverage in t towards the caps, the contribution 
from the horizontal field vanishes. The vertical magnetic field 
strength is related to the integrated thermal electron colunon 
density towards the Galactic poles, DMj_=ne_o^o by 



' 0.812flDMx' 

We note that while there remain controversies on the exact 
values of the scale height Hq and the mid-plane electron den- 
sity nefi, the total column density of thermal electrons is well 
constrained to be ~ 25 pc cm"^ using pulsar DMs at high z 
(see for example, Cordes & Lazio 2003; Gaensler et al. 2008). 

From Eq 13, we found a vertical magnetic field towards 
the north Galactic pole to be consistent with zero (0.00 fiG ± 
0.02 /xG), with a 2>a upper limit on the vertical magnetic field 
strength of 0.07 /xG. On the other hand, the vertical magnetic 
field towards the south Galactic pole is found to be +0.31 /iG 
± 0.03 /xG. We note that the reduced x^s of a model with the 
derived vertical magnetic field strengths towards the Galactic 
poles exceed unity (Xr ~ 7.9 towards the north and Xr ~ 10-6 
towards the south). This is because systematic RM scatter is 
introduced by intrinsic RMs of EGSs and small-scale Milky 
Way foreground electron density and magnetic field fluctua- 
tions. 

Wu et al. (2009) have found, using isothermal magneto- 
hydrodynamic turbulence simulation, a relation between the 
distribution of normalized RM and the line of sight magnetic 
field strength. One can estimate the vertical magnetic field 
strength towards the poles using Equation (3) of Wu et al. 
(2009) and the median and standard deviation of the RM data 
set towards the north and the south caps reported in § 4.5. The 
estimated vertical magnetic field strength towards the north 
and south Galactic pole are roughly jiG and 0.46 /iG re- 
spectively, which agree in general with our results. However, 
as Wu et al. (2009) have pointed out, the relation is valid only 
for a Mach number of unity and thus it is unclear if it holds in 
the diffuse ionized gas at high Galactic latitude. 

There are two pulsars with measured RMs towards the 
Galactic caps: PSRs J0134-2937 and B1237-H25 (Taylor et al. 
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1993; Han et al. 1999). Their coordinates, DMs and RMs are 
listed in Table 3. One can estimate the vertical magnetic field 
between the Sun and the pulsar using DMpuisar and RMpuisar 

g ^ RMpuisar sin 1^7 1 

' 0.812aDMpuisar' 

where b is the Galactic latitude of the pulsar. The vertical 
magnetic field derived from B 1237+25 is +0.044 ± 0.008 
iiG, while that derived from JOl 34-2937 is +0.7 ± 0.1 /iG. 
The vertical magnetic field strength derived using the pul- 
sars roughly agree with results obtained using RMs of EGSs. 
Since PSR J0134-2937 has a much higher DM than PSR 
B 1237+25, it probes through a longer path length in the WIM 
than the northern pulsar 

Taylor et al. (2009) reported a vertical field of -0.14 /iG 
± 0.02 jiG towards the north Galactic pole and +0.30 /uG ± 
0.03 jiG towards the south Galactic pole. While our measured 
vertical field towards the south Galactic cap is consistent with 
Taylor et al. (2009) within errors, our estimation of the verti- 
cal field towards the north Galactic cap disagrees with Taylor 
et al. (2009). This is likely due to the fact that Taylor et al. 
(2009) have averaged over a larger region around the north 
Galactic pole without discarding outliers and anomalous RM 
regions around the pole before performing the fit. Different 
individual RMs derived from the NVSS catalogue and our 
WSRT observations due to multi-RM component sources and 
different ionospheric conditions (§ 4.6) might also contribute 
to this discrepancy. 

5.2.2. Horizontal Field 

In § 5.2.1, we used the fact that the median RM of EGSs 
is zero to reach the conclusion that there is no vertical field 
towards the north Galactic pole. However, it is still possible 
that there exists a horizontal magnetic field at high positive 
Galactic latitude. We can test this by fitting the north Galactic 
cap RMs to a model with only a horizontal component'^ 

RM = -0.812DM_L(BHCOs(£-4)/tan \b\), (15) 

minimizing the ~x} between the observed RMs and the mod- 
eled RMs predicted by the above equation. The best fit pa- 
rameters are Bii=Q-G^) \ l^G and £o= 153°J'[q% with a reduced 
of 7.7. Comparing this model to a model with no mag- 
netic fields using the F-test, the significance of a horizontal 
field of strength 0.6 /uG is only at 2.4 a level. This is not sur- 
prising as no obvious sinusoidal variation of RM as a func- 
tion of I could be seen in the smoothed RM map in the top 
panel of Fig 5. Both the derived horizontal field strength and 
its direction differ from the best fit values towards positive 
mid-latitude (0.39 iiG at 4= 281°) obtained by Taylor et al. 
(2009). As illustrated in Figure 7 of Taylor et al. (2009) , 
the best fit parameters obtained when fitting to high positive 
latitude RMs are different. The authors attributed this to a 
potentially more complicated halo magnetic field. 

We can also test if there is a horizontal magnetic field com- 
ponent in addition to the detected vertical component towards 
the south Galactic pole by performing a least square fit to a 
Galactic magnetic field model with both a vertical and a hor- 
izontal component (Eq 12). The reduced of the best fit 
to such a model is 10.6. The F-test suggests the existence 
of an additional horizontal field to the measured vertical field 
is significant at only 1.7 cr level. This is expected because 

This is obtained by setting B^ = Q jxG in Eq 12. 



no sinusoidal RM variations in t could be seen in the bottom 
panel of Fig 5. We conclude that there is little evidence for a 
horizontal field towards the Galactic poles. 

Instead of attempting to fit for the horizontal field, we can 
subtract its contribution from the measured RMs using the 
best fit values obtained by Taylor et al. (2009) to check if 
it changes our estimation of the vertical magnetic fields in 
§ 5.2.1. This is justified only if the halo magnetic field at 
high latitude is the same as that at mid-latitude. At \b\=ll°, 
the horizontal component contributes a maximum |RM| of ~ 
1.9 rad m"^ towards the north Galactic pole and 4.3 rad m"^ 
towards the south Galactic pole. We have subtracted the con- 
tribution of the Taylor et al. (2009) horizontal halo magnetic 
field from our measured RMs and foimd that it does not al- 
ter the vertical magnetic field estimates presented in § 5.2.1 
within the errors. 

5.3. Galactic Wind 

X-shaped polarization pattern observed in halos of nearby 
edge-on galaxies imply large vertical magnetic fields increas- 
ing with height above the galactic disk (Beck 2008). A kine- 
matic disk dynamo, which generates dominant toroidal mag- 
netic fields, cannot alone explain the existence of these fields. 
As some of these galaxies exhibit evidence of cosmic ray 
driven winds, an alternative explanation of the large vertical 
field is wind advection that transports magnetic fields from 
the disk into the halo, distorting the expected field structure 
from dynamo actions (see for example, Heesen et al. 2009). 

Recent studies by Everett et al. (2008, 2009) have success- 
fully reproduced the Milky Way's diffuse soft X-ray emission 
and synchrotron emission using a ID thermally and cosmic 
ray driven wind model for the Galaxy. In their cosmic ray 
driven wind model, the wind is launched within a Galacto- 
centric radius of ^ 4.5 kpc, but flares to larger Galactocentric 
radii above and below the plane, and hence there is no wind 
launched at the location of the Sun. Similarly, Breitschwerdt 
et al. (1991, 1993) have argued that a wind launched from the 
solar neighborhood would lead to too much cosmic ray es- 
cape and hence inconsistent with the inferred residence time 
of cosmic rays in the Galaxy. Everett et al. (2008, 2009) have 
adopted a flared-cylinder wind geometry where the wind stays 
well confined within a cylinder of constant cross-section up to 
a height Zhreak above/below the Galactic plane. Beyond Zbreak, 
the cross-sectional area of the wind increases as a power law 
of z. At some z, this wind is directly above/below the location 
of the Sun. However, the best fit value of i-break is found to 
be ^ 4 kpc (Everett et al. 2009). At this height above/below 
the Galactic plane, the density of the exponential WIM is too 
low to produce any observable Faraday rotation. Therefore, 
we conclude that the observed RM towards the Galactic caps 
is unlikely to be due to a large scale Galactic wind. 

5.4. Mean Field Dynamo Theory 

The existence of galactic-scale coherent magnetic fields in 
the Milky Way disk and in other normal spiral galaxies with 
significant differential rotation can be explained by the stan- 
dard mean field a-oj dynamo (Beck et al. 1996), although 
the theory and its application to galaxies has been questioned 
on theoretical grounds (see for example Cattaneo & Hughes 
2009, for a recent summary). On a time scale of a few Gyrs, 
this process amplifies and orders the field by turbulence rising 
into the halo, transforming an azimuthal field into a poloidal 
one (the a-effect) and by differential rotation in the underly- 
ing disk, transforming the radial component of the poloidal 
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field back into an azimuthal one (tiie w-effect) (Siiukurov 
2007). In this section, we examine if dynamo theory can pre- 
dict the strength and geometry of a vertical field that is con- 
sistent with the values reported in § 5.2.1. 

The a-u dynamo predicts an azimuthal field that dominates 
over the vertical field by more than a factor of 10 because the 
w-effect operates more efficiently than the a-effect (Ferriere 
2005). We can crudely estimate the expected local vertical 
magnetic field strength by evoking V • B = using a local 
horizontal field of ^ 2 jiG at a pitch angle of 15° (Beck et al. 
1996), and a Galactic disk of total thickness 2h and diameter 
2R (where /z ~ 2 kpc and 7? ~ 15 kpc). The estimated ver- 
tical magnetic field strength at the location of the Sun is ~ 
^radial (h/R) ~ 0.07 /iG. Another estimation of this ratio fol- 
lows from dynamo theory: Bj. ^ Bradiai \/h/R (Ruzmaikin 
et al. 1988). The expected vertical magnetic field strength 
is ~ 0.2 which is similar to the vertical magnetic field 
strength we obtained for the south Galactic cap. The dynamo 
mechanism is thus capable of producing the observed vertical 
magnetic field strength. 

The large scale magnetic field configuration of a galaxy can 
be classified by axial symmetry with respect to the rotation 
axis as well as by vertical synraietry with respect to the galac- 
tic mid-plane. For a strongly differentially rotating galaxy, 
a classical mean-field dynamo favors even field symmetry (a 
quadrupolar field), in which the toroidal component is sym- 
metric across the mid-plane while the vertical component re- 
verses direction (Shukurov 2007). For a weakly differentially 
rotating galaxy, or a halo, odd field synmietry (a dipolar field) 
is preferred, in which the toroidal component reverses direc- 
tion across mid-plane while the vertical component does not 
(Ferriere 2005). However, numerical simulations carried out 
by Ferriere & Schmitt (2000) found comparable growth rates 
of the odd and the even symmetry modes, suggesting that 
the present Galactic magnetic field might be of mixed parity 
rather than of pure even or odd parity. As mentioned in § 1 , 
studies of the Milky Way's field parity have yielded diverg- 
ing results, possibly because they focus on different Galactic 
latitude ranges and because it is difficult to distinguish local 
magneto-ionic effects from a genuine large scale field. The 
most convincing piece of work is the wavelet analysis of all 
sky RMs conducted by Prick et al. (2001). They found that 
the local magnetic field has an even parity - the horizontal 
field component does not reverse direction across the Galactic 
plane. If the Milky Way's large scale field is indeed quadrupo- 
lar in nature, then one expects the vertical magnetic field to 
reverse direction across the mid-plane. On the other hand, 
Han & Qiao (1994) and Han et al. (1997, 1999) concluded 
from an all sky smoothed RM map that the RM distribution 
towards the inner Galaxy is anti- symmetric across the Galac- 
tic plane. The authors attribute this to a large scale dipolar 
field. If this is the case, then one expects the vertical magnetic 
field to have the same direction above and below the Galactic 
plane. In § 5.2.1, we found that the vertical field is consistent 
with zero for b> +1T and that =+0.31 ± 0.03 /xG for b<- 
11° . This is inconsistent with either a quadrupole or a dipole 
large scale field configuration. Even though a simple disk dy- 
namo can produce a vertical field strength comparable to that 
being observed towards the south cap, it cannot explain the 
observed vertical field geometry. 

Sokoloff & Shukurov (1990) demonstrated that the mean 
field dynamo can operate in the halo of a rotating galaxy since 
the mean helicity (the a-effect) of gaseous halos are non-zero 



due to galactic fountains/Parker instabilities. The dominant 
mode is an axisymmetric field with odd parity with respect 
to the Galactic disk; this is of the opposite parity from the 
mode excited by the dynamo action in the Galactic disk. The 
asymmetric RM pattem that we see towards the north/south 
Galactic hemisphere seems to support this idea. In particular, 
the shift of RM structure on the largest scales to h=-\5° found 
by Prick et al. (2001) can be explained by a dipolar field in the 
halo of the Galaxy such that the vertical magnetic field from 
the disk and the halo add up in the southern hemisphere. Sun 
et al. (2008) were able to obtain a reasonable fit to the latitude 
extension of CGPS measurements using an asymmetric halo 
field plus a symmetric disk field. More recently, Taylor et al. 
(2009) found that mid-latitude NVSS RMs can be well fit- 
ted with a ~ 0.4 /xG toroidal halo field that reverses direction 
across the mid-plane. If indeed an anti- symmetric halo field 
and a symmetric disk field co-exist in the Milky Way, then on 
the side of the Galactic disk (in this case ^>0°), the vertical 
component of the halo field and that of the disk field would 
partly cancel out, as they are oppositely directed. These two 
vertical components would add up on the other side of the 
disk (b<Q°). Depending on the relative vertical field strength 
between the halo/disk field and their extension above/below 
the disk, it is not impossible that they can cancel each other 
out exactly towards the northern hemisphere resulting in a net 
RM that is consistent with zero. 

However, this theory has been challenged by numerical 
work. Brandenburg et al. (1992) showed that the halo dynamo 
requires more than a Hubble time to reach the steady state 
configuration and thus the field that one observes at this mo- 
ment might merely be a transient field of mixed parity. Moss 
& Sokoloff (2008) solved the mean field dynamo equations 
for a system with a disk and a halo dynamo. These authors 
were not able to produce any co-existing system of a dipole- 
Uke halo field and a quadrupole-like disk field; instead, one 
always dominates over the other, though Moss & Sokoloff 
(2008) acknowledged that using a larger turbulence diffusivity 
ratio between the halo and the disk might mitigate the prob- 
lem. 

We conclude that the observed vertical magnetic field ge- 
ometry towards the Galactic poles is not consistent with pre- 
dictions from a pure disk dynamo. A separate halo dynamo 
of odd parity could potentially account for the observed verti- 
cal magnetic fields at high latitude but until now no numerical 
simulation has successfully produced a co-existing disk and 
halo fields of opposite parity in a steady state. 

5.5. Primordial origin 

The competing theory to the dynamo origin of galactic 
magnetic fields is that of primordial origin. The primordial 
field theory suggests the following: if the IGM field is frozen 
into the gas, then the total magnetic field would be enhanced 
by a few orders of magnitude as gas clouds collapse to form 
a protogalaxy. This relic field would then be modified by the 
differential rotation of the galaxy, producing the present day 
galactic magnetic fields (see for example. Beck et al. 1996; 
Howard & Kulsrud 1997). The component of the seed field 
parallel to the galactic disk can be removed diffusively and 
through large scale flow, but the component of the magnetic 
field parallel to the rotation axis of the galaxy is trapped (Ruz- 
maikin et al. 1988). 

Since galactic rotation is symmetric with respect to the 
plane, one expects the azimuthal component of such a field 
to reverse its direction across the mid-plane while the vertical 
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component preserves its direction, thus resulting in a dipo- 
lar type field (Ruzmaikin el al. 1988; Beck et al. 1996). The 
strength of a vertical galactic magnetic field of primordial ori- 
gin depends highly on the initial orientation of the intergalac- 
tic seed field with respect to the rotation axis. The collapsed 
gas that forms the proto-galaxy could potentially increase its 
strength by two orders of magnitude (Kulsrud & Zweibel 
2008). This small field is preserved until the present day since 
the total vertical magnetic flux is conserved. The concept of 
a primordial field may be overly simplistic because galactic 
disks probably build up over time through mergers and infall. 
The addition of new material can add new magnetic flux, but 
whatever flux is added is subject to the constraints described 
above. 

The observed vertical field geometry is not consistent with 
that from a pure dipole field of primordial origin: the vertical 
field from the south Galactic pole is directed toward us, while 
the vertical field toward the north Galactic pole is consistent 
with zero. Therefore, we cannot attribute the observed vertical 
magnetic field to a primordial field alone. 

6. TURBULENCE AND RANDOM MAGNETIC FIELD AT HIGH 
GALACTIC LATITUDE 

In § 4.4, we have shown that the RM SFs are flat, imply- 
ing that no specific scale is associated with the observed RM 
pattern towards the north and the south Galactic caps. In this 
section, we will place limits on turbulence properties at high 
1^1 as well as estimate the random magnetic field strength us- 
ing a plane-parallel free electron density model. 

If we assume that the observed RM of an EGS is the sum of 
only its intrinsic RM and the RM produced by the magneto- 
ionic medium in the Milky Way, then the observed structure 
function is the superposition of the intrinsic RM structure 
function and that of the foreground. The SF of intrinsic RMs 
of EGSs is expected to be flat, because intrinsic RM of inde- 
pendent EGSs are uncorrelated with each other, at a saturation 
level of 2(T^jyi i^jji^jjj,. On the other hand, a turbulent Galactic 
foreground with an outer scale l„uier and and an inner scale 
Unner would producc a SF that rises between /inner and /outer and 
saturates beyond lomer at '^o-lMM^gro^mi^ where cxRM.foreground is 
the variance of RM due to random magnetic fields or varia- 
tions in electron density in the foreground ISM. Fig 6 shows 
that the observed RM structure function towards both Galac- 
tic caps saturates at roughly 170 rad^ m""*, which corresponds 
to (TRM,totai ~ 9 rad m . This is then also the upper limit of 

C^RM, foreground* 

The flatness of the observe RM SFs suggests that intrinsic 
RM of EGSs dominate the SF, although one cannot rule out 
the possibility that the outer scale of turbulence in the Galactic 
ISM is below the smallest angular scale probed. This result is 
consistent with previous turbulence studies towards the north 
Galactic pole carried out by Simonetti et al. (1984) and Sun 
& Han (2004) who both found flat RM structure functions on 
scales greater than l°-2° and concluded that the intrinsic EGS 
RM dominates the SF. However, they have reported a satura- 
tion level at log(SFRM) ^ 2.8 which implies a larger total RM 
standard deviation (~ 16 rad m"^) than what we have mea- 
sured. This discrepancy is likely due to the different methods 
used to handle the RM measurement errors in computing the 
SF. It is unclear if the SF of RM measurement errors was sub- 
tracted from the measured SF in these previous works. Also, 
RM outliers in the sample and inaccurate RMs can lead to 
offsets in the SF. We note that the SFs computed by Simonetti 
et al. (1984) and by Sun & Han (2004) were all based on a 



relatively small number of EGSs (52 sources alb > H-60° in 
Simonetti et al. (1984) and 35 sources alh > in Sun & 
Han (2004)), sparsely sampled towards a large region near the 
pole. 

Next, we use the coherent vertical magnetic field strength 
derived in § 5.2.1 and properties of turbulence in the Galac- 
tic halo to estimate the random magnetic field strength at 
high Galactic latitude. We construct a cell model similar to 
those used to estimate the random field in the inner Galaxy 
(Gaensler et al. 2001), the Large Magellanic Cloud (Gaensler 
et al. 2005) and flie Small Magellanic Cloud (Mao et al. 2008). 

We assume that the WIM is composed of plane -parallel 
cells of size /outer in pc, equivalent to the outer scale of tur- 
bulence. We also assume that all the sight lines towards EGSs 
in a 13° cone around the Galactic poles are parallel. We in- 
troduce a clumpy WIM, where the volume average electron 
density ngiz) (Eq 1 1) is the product of the filling factor (the 
fraction of the line of sight at height z filled with electrons) 
and the internal cell electron density N{z). We have adopted 
the evolution of the filling factor and the ceU electron density 
with height above the Galactic plane by Gaensler et al. (2008). 

/(^) = /oe^l'l/''^ (16) 

N{z) = N(>e-\'^l"\ (17) 

where /o=0.04 and A'b=0.34 are the mid-plane filling factor 
and internal cell density; Hj ^ 700 pc and Hf^ ^ 510 pc are 
the scale heights of the filling factor and the internal cell den- 
sity respectively. Eq 16 and 17 were derived by Gaensler 
et al. (2008) using independent constrains on tigiz) from pul- 
sar DMs and on rt^iz) from WHAM Ha measurements to- 
wards the Perseus arm. For z > 1 .4 kpc, beyond the height 
of WIM observed by WHAM towards the Perseus arm, the 
above equations are no longer valid as n^(z) is not constrained. 
Below z of 1.4 kpc, the filUng factor increases exponentially 
with height as one expects the WIM to dominate over the neu- 
tral component at higher z- In each cell, we assume that there 
is a constant coherent magnetic field B,sinb k, (for b ss 90°) 
parallel to the line of sight and a magnetic field with constant 
strength Bran but random orientations with respect to the Une 
of sight in different cells. For simplicity, we assume that 
and Bran are independent of z. The line of sight magnetic field 
through a particular cell m is given by 

cos6. (18) 

The average RM across many random sight lines through M 
cells is given by 

{RMMceiis)=0.Sl2B, f neiz)dz, (19) 

where M * /outer=l-4 kpc. On the other hand, the average of 
squared RM over many random sight lines depends explicitly 
on the cell electron density and the filling factor. 
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(20) 

Using Equation 19, 20 and the definition of standard devi- 
ation, one can write down the expression of the random mag- 
netic field strength in the WIM in terms of the standard devi- 
ation of RM arised in the foreground magneto-ionic medium. 
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CRM, foreground 



Br. 



(21) 

We have placed an upper limit to the the standard devia- 
tion of RM due to the Galactic foreground CTRMjoreground of 
approximately 9 rad m"^ earlier in this section. We assume 
an outer scale of turbulence (/outer) of ~ 90 pc estimated us- 
ing the autocorrelation function of synchrotron radiation mea- 
sured towards the north Galactic pole (Lazaryan & Shutenkov 
1990). As Eq 16 and 17 are physically meaningful only for 
z<1.4 kpc, we choose to carry out the integral in Eq 21 to z = 
M * /outer =1.4 kpc. Using a uniform magnetic field strength 

= (UG, we find that Brm < 1 -5 jjG towards the north Galac- 
tic pole. Assuming that turbulence towards the south Galactic 
pole has the same outer scale as that towards the north and 
using a uniform vertical field of \B.\ = 0.31 /iG, we obtain an 
upper hmit of the random magnetic field strength of 1.4 jiG 
towards the south Galactic pole. Our estimated random field 
strength of ^ 1 is systematically smaller than other es- 
timates of the random field in the Milky Way disk: 4-6 fiG 
was inferred from DM and RM of pulsars (Ohno & Shibata 
1993) and 5 ± 2 /iG inferred using the Haslam 408 MHz syn- 
chrotron intensity map (Beck et al. 1996; Shukurov 2007). If 
one assumes energy equipartition between the turbulence and 
random magnetic fields, Bran is approximately 

Bran = , (22) 

where v is the rms turbulent velocity and p is the density (Beck 
et al. 1996). Since both the density and the turbulent velocity 
in the halo is smaller than that in the Galactic disk, the random 
magnetic field in the halo is expected to be smaller than the 
values estimated for the Galactic plane. Cox (2005) assumed 
equipartition between thermal and non-thermal pressure and 
derived a field strenth of ^ 2 ^G at a few kpc above/ below 
the Galactic plane. This is consistent with the ^ l/xG random 
magnetic field that we computed. 

7. SUMMARY AND OUTLOOK 

In this paper, we have presented an RM survey with the 
ATCA and the WSRT of polarized extragalactic radio sources 
towards the north and the south Galactic poles at \b\ > 11° . 
Using rotation measure synthesis, we have obtained 813 reli- 
able RMs towards the Galactic poles. No preferred RM fiuc- 
tuation scale was apparent from the flat RM structure func- 
tions towards the Galactic caps. After discarding outliers and 
anomalous RM regions in § 4, we obtain a median RM of 0.0 
± 0.5 rad m"^ towards the north Galactic pole; and a median 
RM of H-6.3 ± 0.7 rad m"^ towards the south Galactic pole. 

In § 5.1, we have ruled out the possibiUty that local 
sources/events such as the LISM, the LB and a Parker's in- 
stability loop produce the observed RM pattern. In § 5.2.1 
and § 5.2.2, we have derived the halo magnetic field properties 
from the observed RMs assuming that they are produced by 
the diffuse interstellar free electrons. We found no evidence 
for vertical and horizontal magnetic field towards the north 
Galactic pole. On the other hand, a vertical field of strength 
H-0.31 ± 0.03 /iG was detected at >9 a towards the south 
Galactic pole, but there is no evidence for an additional hor- 
izontal component. Although a dynamo or a primordial field 
can explain the derived vertical magnetic field strength to- 
wards the Galactic poles, a pure dipole/ quadrupole field can- 
not explain the geometry of the observed vertical field across 



the mid-plane. One possible explanation of the derived mag- 
netic field properties is that proposed by Sokoloff & Shukurov 
(1990), in which a disk and a halo dynamo of different parities 
are simultaneously at work in the Galaxy. This could poten- 
tially lead to part cancellation of RM produced by the verti- 
cal magnetic field in the northern Galactic hemisphere, which 
is compatible with our RM measurements. However, until 
now, no numerical simulation has been able to produce a co- 
existing system of a dipole-hke halo field and a quadrupole- 
hke disk field. Numerical works that explore larger parameter 
space (especially turbulent diffusivity) is needed to test this 
hypothesis. Finally, we have estimated the random magnetic 
field strength in the halo of the Milky Way by constructing a 
plane-parallel cell model of the WIM and the standard devia- 
tion of RMs to derive a random magnetic field strength of 
1 /iG in the Galactic halo, which is smaller than the random 
field in the mid-plane of the Galaxy, but in equipartition with 
the lower turbulent energy density inferred for the halo. 

Exploration of cosmic magnetism is one of the key sciences 
of the next generation radio telescopes - the Square Kilometre 
Array (SKA) and its prototypes such as the Australian Square 
Kilometre Array Pathfinder (ASKAP), which are capable of 
providing accurate RMs of EGSs densely sampled over the 
entire sky (Johnston et al. 2007, 2008). For example, one of 
the approved Survey Science project of ASKAP - the Polar- 
ization Sky Survey of the Universe's Magnetism (POSSUM) 
aims to perform RM synthesis and obtain a grid of RMs over a 
large fraction of the sky. Similar projects will provide a more 
detailed picture on the magnetic field and turbulence proper- 
ties at high Galactic latitude. Also, pulsar searches at high 
Galactic latitudes will allow one to probe the vertical mag- 
netic field as a function of height above/ below the Galactic 
disk by simultaneously using pulsar DM and RM, which can 
further constrain the structure of the Galactic halo magnetic 
field. 
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Least square tit EM = -3.0 +/- 1.2 rad/m^ RM Synthesis RM = -3.1 +/- 1.2 rad/m^ 

//■/ = 1.91 SNR = 141.3 
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Fig. 1 .— The least square fit and RM synthesis results for an EGS towards the north (RA=12:43:41.17, DEC=+31:24:17.88) in the top row and an EGS towards 
the south (RA=00:31:02.73, DEC=-22:07:06.50) in the bottom row. The left panel shows the position angle of linear polarization versus relation for each 
source, with the best fit plotted as a dotted line. The right panel is the amplitude of the linear polarization as a function of Faraday depth, the thin dotted line is 
the dirty spectrum and the thick solid line is the deconvolved spectrum. 
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Fig. 2. — Histogram of the RM distribution towards tlie north Galactic pole in the top panel and the south Galactic pole in the bottom panel, binned every 
4 rad m"^. The histograms include all raw RMs (solid outlined histograms) as well as the final RMs after removing extrema and anomalous regions (shaded 
histograms). The thick vertical lines indicate the 3.2 cr cut-off boundary for the extreme RM rejection scheme. The arrow in the top panel indicates the median 
(0.0 rad m"^) of the final north Galactic cap RM data set, while the arrow in the bottom panel indicates the median (+6.3 rad m"^) of the final south Galactic cap 
RM data set. The standard deviations of the final data sets are 9.2 rad and 8.8 rad m"^ towards the north and the south Galactic pole respectively. 
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Fig. 3. — Rotation measure distribution overlaid on tlie Ha emission map (Finlcbeiner 2003) towards the north Galactic pole in the top panel and towards the 
south Galactic pole in the bottom panel (after filtering sources with low signal-to-noise ratio or high polarized fractions; sources with extreme value of RMs or in 
anomalous RM regions are still shown). The color scale is in units of Rayleighs. Positive (Negative) RMs are denoted by filled (open) circles with their diameters 
proportional to the magnitude of RM. The largest circle corresponds to a RM value of +93 (+76) rad m"^ towards the north (south) Galactic pole. Sources with 
RMs consistent with zero at 1 cr are denoted by asterisks. The blue cross in the top panel indicates the centre of the Coma cluster. The dotted green curve in the 
bottom panel marks the projection of the LB wall: sight lines to the east of the boundary penetrate the LB wall. The resolution change at the bottom of the top 
image corresponds to the change of Ha data being used to create the composite map from the WHAM to the Southern H-Alpha Sky Survey Atlas (SHASSA), 
the latter has better spatial resolution but much poorer sensitivity. The Ha image towards the south Galactic pole is a combination of SHASSA and WHAM data 
in most of Galactic Quadrant 2 and 3 and it consists of SHASSA data only in Quadrant 1 and 4. 
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Fig. 4. — The scatter plot of EMs (from WHAM) against RMs along sight lines to EGSs towards the north (top panel) and the south (bottom panel) Galactic 
caps. All EOS sight lines towards the north Galactic pole have corresponding WHAM EM measurements, while only 192 sight lines towards the south Galactic 
pole has corresponding WHAM EM measurements. For clarity, measurement errors associated with EM and RM are not plotted. No clear trend of correlations 
between EM and RM can be seen. 
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Fig. 5. — The smooth RM map towards the noi1h (top panel) and the south (bottom panel) Galactic cap in SIN projection. The Galactic poles are located at 
the center of the images. The median RMs in 2° by 2° pixels are computed. Blanked regions represent pixels with insufficient number of EGSs (<3 ) to obtain 
reliable statistics. The color scale to the right of the figure is in units of rad nC^ and is chosen such that rad m"- is represented by white. Pixels values that 
deviate more than 1 .64 cr from the median of the entire RM distribution are regarded as anomalous regions and are indicated by green rectangles. One such cell is 
identified towards the north Galactic cap and consequently, the 3 RMs within the cell are discarded for the puipose of studying the large scale Galactic magnetic 
field. No anomalous RM region is identified towards the south Galactic cap. 
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Fig. 6. — The RM structure function towards the north (top) and the south (bottom) Galactic cap computed using Eq 8 and 9. The structure function has been 
binned in equal log interval (0.25) and the minimum number of correlations per bin is 10. The error bars denote the standard error in the mean of each bin. 
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South Galactic Pole 




Fig. 7.— The comparison between NVSS RMs derived by Taylor et al. (2009) and our RMs. NVSS RMs are plotted against WSRT (ATCA) RMs towards 
the north (south) Galactic pole in the upper (lower) panel. The solid line indicates where NVSS RMs and our RMs agree with each other. For clarity, error bars 
associated with the measurements are not drawn. 
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TABLE 3 

Radio Pulsars with measured RMs at \b\ > 11° 



Name 




b 


DM (pc cm-^)" 


RM(rad m-2) 


B 1237+25 


252.45 


86.54 


9.242 ± 0.006 


-0.33 ± 0.06 


J0134-2937 


230.25 


-80.25 


21.806 ± 0.006 


13±2<= 



" Hobbs et al. (2004) 
Taylor etal. (1993) 
Han et al. (1999) 



